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ly-^ ■ Abstract. Two-body bound states such as mesons are described by solutions of the Bethe- 

Salpeter equation. We discuss recent results for the pseudoscalar and vector meson masses and 
leptonic decay constants, ranging from pions up to cc bound states. Our results are in good 
agreement with data. Essential in these calculation is a momentum-dependent quark mass 
^ ' function, which evolves from a constituent-quark mass in the infrared region to a current- 

quark mass in the perturbative region. In addition to the mass spectrum, we review the 
electromagnetic form factors of the light mesons. Electromagnetic current conservation is 
manifest and the influence of intermediate vector mesons is incorporated self-consistently. The 
results for the pion form factor are in excellent agreement with experiment. 



1. Dyson Schwinger equations 

O \ The set of Dyson-Schwinger equations form a Poincare covariant framework within which to 

study hadrons [1, 2]. In rainbow-ladder truncation, they have been successfully applied to 
calculate a range of properties of the light pseudoscalar and vector mesons, see Ref. [2] and 
references therein. 

^ ' The DSE for the renormalized quark propagator S{p) in Euclidean space is [1] 

s{p)-^ = iZ2{C)i> + z^{0m{0 + Zi(c) j-^a^^Ap - 1) h^^S{q)ri{q,p) , (1) 

where Dfj,^{p — q) and T1{q;p) are the renormalized dressed gluon propagator and 
quark-gluon vertex, respectively. The most general solution of Eq. has the 

form S(p)^^ = i^A(p'^) -\- B{p'^), renormalized at spacelike according to A{('^) = 1 and 
B{('^) = m{C,) with m{C,) the current quark mass. 

Mesons are described by solutions of the homogeneous Bethe-Salpeter equation (BSE) 

Th{p+,P-;P)= J-^K{p,q;P)S{q+)TH{q+,q^;P)S{q.), (2) 

at discrete values of = —Mfj, where Mh is the meson mass. In this equation, pj^ = p + P/2 
and p^ = p — P/2 are the outgoing and incoming quark momenta respectively, and similarly for 
q±. The kernel K is the renormalized, amputated qq scattering kernel that is irreducible with 
respect to a pair of qq lines. Together with the canonical normalization condition for qq bound 
states, Eq. Q completely determines the bound state Bethe-Salpeter amplitude (BSA) Th- 



Different types of mesons, such as pseudoscalar or vector mesons, are characterized by different 
Dirac structm'es. 

To solve the BSE, we use the ladder truncation, 

K{p, q- P) ^ -47ra^f^((^, - qf) D^^ip - q)^j^ ^ ^7, , (3) 

in conjunction with the rainbow truncation for the quark DSE, Eq. T]^{q,p) ^^Iv ^-^d 
Zig'^D^j,y{k) A-KQ.'^^ {k'^)D^^^,{k) , with k = p — q. Here, D^^{k) is the free gluon propagator 
in Landau gauge, and a'^^{k'^) an effective quark-quark interaction, which reduces to the one- 
loop running coupling of perturbative QCD (pQCD) in the perturbative region. This truncation 
preserves both the vector Ward-Takahashi identity (WTI) for the qq'j vertex and the axial- vector 
WTI, independent of the details of the effective interaction. The latter ensures the existence 
of massless pseudoscalar mesons associated with dynamical chiral symmetry breaking (D^SB) 
in the chiral limit [3,4]. In combination with an impulse approximation, the former ensures 
electromagnetic current conservation [5,6]. 



2. Quark propagator: an overview 

A momentum-dependent quark mass function M(p^) = Bijp') /Aijp') is central to QCD. In the 
perturbative region this mass function gives the one-loop perturbative running quark mass 
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with the anomalous mass dimension 7.^ = 12/(33 — 2A^j). Dynamical chiral symmetry breaking 
means that this mass function is nonzero even though the current-quark masses are zero. In the 
chiral limit the mass function is [7] 
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with {qqY the renormalization-point-independent vacuum quark condensate [3]. 

It is a longstanding prediction of DSE studies in QCD that the dressed quark propagator 
receives strong momentum-dependent corrections at infrared momenta, see e.g. Refs. [1, 2] 
and references therein. Provided that the (effective) quark-quark interaction reduces to the 
perturbative running coupling in the ultraviolet region, it is also straightforward to reproduce 
the asymptotic behavior of Eqs. (0J) and © [8,9]. Both these phenomena are illustrated in the 
left panel of Fig. ^ In this figure one can also see that the dynamical mass function of the u and 
d quarks becomes very similar to that of quarks in the chiral limit in the infrared region. This 
is a direct consequence of D^SB, and leads to a mass function of the order of several hundred 
MeV for the light quarks in the infrared, providing one with a constituent mass for quarks inside 
hadrons, even though the corresponding current quark masses are only a few MeV. 

These predictions were recently confirmed in lattice simulations of QCD [10-12]. Quantitative 
agreement between the lattice simulations and the DSE results for the quark propagator 
functions can be obtained within the rainbow truncation via a suitable choice for the effective 
quark-quark interaction [13]. Pointwise agreement for a range of quark masses requires this 
interaction to be flavor-dependent [14], suggesting that dressing the quark-gluon vertex r^(g,p) 
is important. Indeed, both lattice simulations [15] and DSE studies [16-18] of this vertex indicate 
that T\,{q,p) deviates significantly from a bare vertex in the nonperturbative region. A (flavor- 
dependent) nonperturbative vertex dressing could make a significant difference for the solution 
of the quark DSE [16, 19], as can be seen from the right panel of Fig. ^ The consequences of a 
dressed vertex for the meson BSEs are currently being explored [16] and indications are that in 
the pseudoscalar and vector channels, the effects are small [16,20,21]. 
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Figure 1. The dynamical quark mass function MijP') for different quark flavors (left, 
adapted from Ref . [3] ) and a comparison of the chiral-limit quark mass function with the lattice 
data [11,12] (right, adapted from Ref. [19]). 



3. Pseudoscalar mesons: ground and excited states 

The meson spectrum contains three pseudoscalars with quantum numbers I^{J^)L = 1^(0^)5 
and masses below 2GeV: 7r(140); 7r(1300); and 7r(1800). In a constituent-quark model, these 
mesons are viewed as the first three members oi a qq n ^Sq trajectory, where n is the principal 
quantum number with the ground state ttq (the 7r(140)), and the others are its first two radial 
excitations, vri and 7r2. The pseudoscalar trajectory is particularly interesting, because its lowest 
mass member is QCD's Goldstone mode. Therefore an explanation should describe both chiral 
symmetry and its dynamical breaking as well as a correlation of the ground and excited states 
via an approximately linear radial Regge trajectory. The latter is easily realized in Poincare 
invariant quantum mechanics [22] but the former is not. 

3.1. Chiral symmetry 

The chiral properties of QCD are manifest in the axial-vector WTI, which reads 

P^r5^(g+,g_;P) = S-\q+)n5+n5S-\q-)- 2img{C)'r5iq+,q-;P), (6) 

with r5^((7+, g_; P) and T^{q^, q^; P) the renormalized dressed axial- vector and pseudoscalar 
vertices, each satisfying an inhomogeneous extension of Eq. Equation © is an exact 

statement in QCD implying that the kernels of the quark DSE, Eq. (^, and of the BSEs have 
to be intimately related. A weak coupling expansion of the DSEs yields perturbation theory 
and satisfies this constraint, but is not useful in the study of intrinsically nonperturbative 
phenomena. However, a systematic and symmetry preserving nonperturbative truncation 
scheme exists [16, 20, 21], allowing for both elucidation and illustration of the consequences 
of the axial- vector WTI. 

Pseudoscalar mesons appear as pole contributions to the axial-vector and pseudoscalar 
vertices at = —Ml . The residues of these poles are 

f.^P, = Z2(C)/^Tr[757M5(g+)r^„(<7+,g_;P)S(g_)] , (7) 
ip^SC) = ^4(0 J^TT[^,S{q+)Kjq+,q.;P)S{q.)], (8) 



both of which are gauge invariant and cutoff independent. It follows from Eq. ® that these 
residues satisfy the following exact identity in QCD [4, 23] 



U„Ml^ = 2m(C)p.„(C), (9) 

valid for every 0~ meson [24], irrespective of the magnitude of the current quark mass [25]. 

For the ttq, DxSB yields: limm^o /ttq / and limm_,o piC) = -'^{m)^ fno 0. Hence, 
the Gell-Mann-Oakes-Renner relation emerges as a corollary of Eq. Q and the ground state 
pion is QCD's Goldstone mode [4]. For the n > 1 pseudoscalar mesons one has M^r^^^ > M^g 
by assumption, and hence M,r„>o 7^ in the chiral limit. Furthermore, the ultraviolet behavior 
of the quark-antiquark scattering kernel in QCD guarantees that PwniC) is finite in the chiral 
limit. Hence, it is a necessary consequence of chiral symmetry and the axial-vector WTI that 
/tt^ vanishes in the chiral limit for all excited pseudoscalar mesons (i. e. for all n > 0) [24]. 

3.2. Numerical results 

We now illustrate the exact results reviewed above in a model that both preserves QCD's 
ultraviolet properties and exhibits D^SB, namely the rainbow-ladder truncation of the set of 
DSEs. For the infrared behavior of the effective quark-quark interaction, a^^{k'^), we employ 
an Ansatz [3, 26] that is sufficiently enhanced in the infrared to produce a realistic value for the 
vacuum quark condensate of about (240 GeV)^. The model parameters, along with the quark 
masses [26,27], are fitted to give a good description of the chiral condensate, M^jj^j^^ and f-,^. 
The obtained quark propagator functions agree qualitatively with lattice simulations. 

The main results for the pseudoscalar and vector meson masses and leptonic decay constants 
are summarized in Tableland illustrated in Fig. El Regarding the table one should note that the 
rainbow-ladder truncation gives "ideal" flavor mixing and the study uses rhu = w-d 7^ rhg 7^ m-c- 
Thus the first row simultaneously describes four degenerate mesons for each column, namely, 
a {tid, uu — dd,du} isotriplet and a uu + dd isosinglet. The second row describes ss mesons, 
and the third row cc mesons. For vector mesons ideal mixing is a very good approximation, 
though this is not the case for the pseudoscalar ground states. However, the experimental 
degeneracy of the 7r(1300) and ry(1295) suggests [28] that ideal mixing is almost realized for the 
excited pseudoscalar mesons. This supports the interpretation of the ry(1295) and r/(1470) as 
the radial excitations of the 77(548) and r7'(958), with quark content almost entirely uu + dd and 
ss respectively [27]. Figure 121 shows bound-state masses and leptonic decay constants for the 
ground and excited states of pseudoscalar mesons as well as ground-state vector-meson masses 
as functions of the current-quark mass. It illustrates that 772,7^0 vanishes in the chiral limit, while 
rrip and m.,^-^ do not; on the other hand, /t^^ vanishes in the chiral limit, while /ttq does not. 



Table 1. Results for the qq 0~ ground and first radially excited states as well as qq 1^ ground 
states. The current-quark masses are = 5.4 MeV, nig = 124MeV, and rUc = 1.34 GeV at 

C = 1 GeV, and all quantities are given in GeV; for the leptonic decay constants we follow the 
conventions of Ref. [25,26]. Experimental data are taken from [28] except for /^^ [29]. 
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Figure 2. The masses M of the pseudoscalar ground and first radially excited states (left) and 
the corresponding weak decay constants (right) as functions of the current-quark mass niq. For 
comparison, we also include the mass of the ground-state vector mesons. 

Satisfaction of the axial- vector WTI can be checked by virtue of Eq. Q for any pseudoscalar 
meson in Table ^ For the ground and excited states the inaccuracies are below 1% and 5%, 
respectively. The latter are larger because of: (1) the need to project out the ground state in 
order to calculate the excited state, (2) the smallness of close to the chiral limit, and (3) the 
larger bound-state mass, leading to a larger region of analytic continuation for the quark DSE 
solution in the complex plane. 

4. Electromagnetic form factors 

Meson-meson-photon interactions can be described in generalized impulse approximation by 

r''^(P,Q,i^) = |-|^Tr[5"(g)^"^(<7,<7^P)5^gO^''^"(9^9'^Q)5^(9'0^^'(9^9;i^)] , (10) 

where q — q' = P,q' — q" = Q,q" — q = K, and momentum conservation dictates P + Q + K = 0. 
In Eq. ()1U() . 5* is the dressed quark propagator with flavor index i, and T^^{k, k'; P) stands for 
a generic vertex function with incoming quark flavor j and momentum k', and outgoing quark 
flavor i and momentum k. Depending on the specific process under consideration, this vertex 
function could be a meson BSA or a quark-photon vertex. 

The quark-photon vertex, r^(p+,p_; Q), with Q the photon momentum and p± the quark 
momenta, is the solution of the inhomogeneous BSE 

T^{p+,p_;Q) = Z2{0l>. + J-^Kip,q;Q) S{q+)T^{q+,q_;Q) S{q.) . (11) 

Solutions of the homogeneous version of Eq. (|llj) define vector meson bound states at timelike 
photon momenta = —My. It follows that T^{p^,p^; Q) has poles at these locations [6] (see 
also the vector-meson masses in Table 

4.1- Pion and kaon elastic form factors 

There are two diagrams that contribute to meson electromagnetic form factors: one with the 
photon coupled to the quark and one with the photon coupled to the antiquark respectively. 
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Figure 3. Left: The pion form factor Q'^F{Q'^), compared to the experimental results from 
Refs. [30,31]. Right: The transition form factor F^^^{Q'^) together with experimental data [40]. 

With photon momentum Q, and incoming and outgoing meson momenta P^Q/2, we can define 
a form factor for each of these diagrams [6] 

2P,F,5^5(Q2) = lf~\P-Q/2,Q,-{P + Q/2)). (12) 

We work in the isospin symmetric limit, and thus F.,^{Q'^) = Fuu-uiQ"^)- The and form 
factors are given by Fk+ = |-F«s;m + \Fus-s and F^o = -\Fds-d + \Fds;s, respectively. 

Our result for Q'^Ft^ are shown in the left panel of Fig. |31 In the timelike region, and in 
the spacelike region up to about Q2 = 2GeV^ F^ can be described very well by a monopole 
with our calculated p-mass, nip = 742 IVIeV (note that our calculated p-mass is slightly below 
the experimental value). Above this value, our curve starts to deviate more and more from this 
naive vector meson dominance (VMD) monopole. Our result is in excellent agreement with the 
most recent JLab data [30]; it would be very interesting to compare with future JLab data in 
the 3 to 5 GeV^ range, where we expect to see a significant deviation from the naive monopole 
behavior. The pQCD prediction is significantly smaller than our results at 4 GeV^; we anticipate 
that true perturbative behavior sets in somewhere between 10 and 20 GeV^ [32]. 

Recent lattice simulation also indicate that in the region between and 2 GeV^ the pion 
form factor can be represented by a VMD-like monopole [33,34]. This VMD-like behavior of 
the pion form factor in the spacelike region appears to be valid almost independent of the pion 
mass, both in our calculations and in the lattice simulations, though a monopole fit results in a 
VMD-mass which is slightly less than the actual vector meson mass. Current lattice simulations 
are not accurate enough, nor do they extend to large enough values of Q^, to detect a significant 
deviation from VMD-like behavior above 2 GeV'^. 

Also our results [6] for Fx agree quite well with the available experimental data, as do both 
the neutral and the charged kaon charge radius [35,36], see Table El It should be noted here 
that the form factor is obtained by taking the difference of two numbers, Fds-diQ"^) and 
Fds-siQ"^), which are both close to one for near zero. It is therefore much more sensitive to 
details of the model: both kaon loops and a flavor dependence of the (effective) quark-quark 
interaction will have a significantly bigger effect on the form factor than on the form 
factor. This could explain why it does not agree as well with experiment as our other results. In 
this respect one should note that the absolute deviation between our calculated charge radius 
and the experimental charge radius is similar for the charged and neutral kaon. 



4-2. Radiative vector-meson decays 

We can describe the radiative decay of the vector mesons using the same loop integral, Eq. (fTU]) . 
this time with one vector meson BSA, one pseudoscalar BSA, and one 557- vertex [37]. The 
generic structure of a vector-pseudoscalar-photon vertex is 

I-l^-(P,Q,-{P + Q)) = e^,o.pPaQpfab;aiQ^), (13) 

where P is the vector momentum and Q the photon momentum. The on-shell value gives us the 
coupling constant, and can be used to calculate the partial decay width of the vector mesons. 
For virtual photons, we can define a transition form factor Fyp^(Q^), normalized to 1 at = 0, 
which can be used in estimating meson-exchange contributions to hadronic processes [38,39]. 
In the isospin limit, both the ir^ 7 and vr^ 7 vertices are given by 

^e^^^pPaQpF^^iQ^) = ll-^'-{P,Q,-{P + Q)). (14) 
mp 

The ujTT'j vertex is a factor of 3 larger, due to the difference in isospin factors; however, the form 
factor Fyp^{Q'^) is the same for pir^ and ujtt^. In contrast to the elastic form factors, this 
transition form factor falls off significantly faster than a VMD-like monopole, as can be seen in 
the right panel of Fig. |21 Only in the timelike region, near the vector meson pole, do we see a 
true VMD-like behavior. 

As Eq. ()14() shows, it is gvP-y/iTT-v that is the natural outcome of our calculations. Therefore, 
it is this combination that we give in Table [2 together with the corresponding partial decay 
widths [37]. As anticipated, the partial decay width u; — s- 7r7 is indeed (approximately) nine 
times larger than the u tt'j partial decay width. Note that part of the difference between 
the experimental and calculated decay width comes from the phase space factor because our 
calculated vector meson masses deviate up to 5% from the physical masses. 

For the K*K^ decays (and corresponding form factors) [41] the situation is more complicated 
owing to the interference of the diagrams with the photon coupled to the s-quark and to the 
u- or d-quark. In the SU(3) flavor limit, the charged K* K 'j vertex becomes equal to the pTT'j 
vertex, whereas the neutral K* K 7 vertex is twice as large in magnitude 

i^-^^X n*K,{Q^) = lUAQ') - lUAQ') , (15) 

V mx* J 

^K^K.iQ') = -\h-sAQ') - yu-sAQ') ■ (16) 

In Table 121 we see indeed that the partial decay width of the neutral K* — s- is larger than 
that of the charged K* — > X7, though not by a factor of four as it would in the SU(3) flavor 



Table 2. Overview of our results for the pseudoscalar-meson charge radii squared, all in fm^, 
with an estimated combined numerical error of less than 0.01 fm^. Also included are our results 
for the vector meson radiative decays: Ty^p^ in keV and g/m in GeV~^. 







Tp± g/m 


F^ g/m 


F^*± g/m 


F^*o g/m 


calc. 
expt. 


0.44 0.38 -0.085 
0.44 0.34(5) -0.052(26) 


53 0.69 
68 0.74 


479 2.07 
757 2.38 


90 0.99 
50 0.84 


130 1.19 
116 1.27 



limit. Furthermore, we see that the deviation between our calculation and experiment is largest 
for the charged K* — > decay. Again, this can be understood since this decay is sensitive 
to the difference between the impulse diagrams, and therefore more sensitive to details of the 
model and its omissions. 
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